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COVER DNA microarrays for analyzing complex DNA samples. Shown is a two-color fluorescent scan of 
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GENOME METHODS N °p c * This meter*' n» y e-» 

A DNA Microarray System for Analyzing 
Complex DNA Samples Using Two-color 
Fluorescent Probe Hybridization 

Dari Shalon, 1 ' 4 Stephen J. Smith/ and Patrick O. Brown 1 - 2 - 5 

toward Hughes Medical Institute and Departments of biochemistry and ^Molecular and Cellular 

Physiology, Stanford University, Stanford, California 94305 



Detecting and determining the relative abundance of diverse individual sequences in complex DNA samples is 
a recurring experimental challenge in analyzing genomes. We describe a general experimental approach to 
this problem, using microscopic arrays of DNA fragments on glass substrates for differential hybridization 
analysis of fluorescently labeled DNA samples. To test the system, 864 physically mapped X clones of yeast 
genomic DNA, together representing >7S% of the yeast genome, were arranged into 1.8-cm x 1.8<m arrays, 
each containing a total of 1744 elements. The microarrays were characterized by simultaneous hybridization 
of two different sets of isolated yeast chromosomes labeled with two different fluorophores. A laser 
fluorescent scanner was used to detect the hybridization signals from the two fluorophores. The results 
demonstrate the utility of DNA microarrays in the analysis of complex DNA samples. This system should 
find numerous applications in genome-wide genetic mapping, physical mapping, and gene expression studies. 



Many problems in genome analysis depend on 
determining what specific sequences are repre- 
sented in a complex DNA or RNA sample and at 
what abundance, for example, what genes are 
represented in a specific chromosome band or 
YAC clone, what intervals are amplified or de- 
leted in a particular cancer cell or what genes are 
expressed in specific cells under specific condi- 
tions. As a general approach to this problem, we 
have developed a system for making microarrays 
of DNA samples on glass substrates, probing 
them by hybridization with complex fluorescent- 
labeled probes, and using a laser-scanning micro- 
scope to detect the fluorescent signals represent- 
ing hybridization. Fluorescent labeling allows for 
simultaneous hybridization and separate detec- 
tion of the hybridization signal from two or more 
probes. This in turn allows very accurate and re- 
liable measurement of the relative abundance of 
specific sequences in two complex samples. 

RESULTS 

Array Hybridization Pattern 

Figure 1 shows the two-color fluorescent scan of 
a yeast genomic array following hybridization 

*»r**ent address: Syntenl. Inc. Palo Aho, California 94305. 
'Corresponding author. 

E-MAIL p brown cm9m.fttanf0rd.adu, http://cmgm. 
s tanf ord.edu/pbrown; FAX (4 1 S) 72 J* 1 399. 



with a mixed probe consisting of lissamine- 
labeled DNA from the 6 largest yeast chromo- 
somes together with fluorescein -labeled DNA 
from the 10 smallest yeast chromosomes. A red 
color indicates that yeast sequences present in 
the lissamine-labeled hybridization probe hy- 
bridized to an array element. A yellow-green 
color indicates that yeast sequences present in 
the fluorescein-labeled hybridization probe hy- 
, bridized to an array element. An orange color in- 
dicates cross-hybridization of both chromosome 
pools to an. array element (e.g., dispersed repeti- 
tive elements, such as Tyl elements). 

Each clone was spotted twice, resulting in du- 
plicate hybridization patterns in adjacent quad- 
rants of the array. Control DNA spots, which 
were randomly amplified in the same manner as 
the X clone array elements, are located in the bot- 
tom corner of each quadrant. "A" points to a pair 
of spots containing total yeast genomic DNA. 
These spots appear orange because both chromo- 
some pools hybridized to yeast genomic DNA. 
The negative controls are as follows: "B" points 
to a pair of spots of wild-type X DNA, "C" points 
to a pair of human genomic DNA spots, and "D" 
points to a pair of 6X174 DNA spots. The lack of 
a hybridization signal at these three negative 
control spots indicates that the hybridization was 
specific for yeast sequences. 
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Figure 1 Two-color fluorescent scan of a 1 .8-cm x 1 .8-cm yeast array 
of \ clones of yeast genomic DNA. The DNA spots are spaced "at a 
distance of 380 jim from center to center. A probe mixture consisting of 
DNA from the 6 largest yeast chromosomes (4, 7, 12, 1 3, 1 5, 1 6) labeled 
with lissamine (red dots) and DNA from the 1 0 smallest yeast chromo- 
somes (1, 2, 3, 5, 6, 8, 9, 10, 11, 14) labeled with fluorescein (yellow- 
green dots) was hybridized to the array. A pair of yeast genomic DNA 
spots (A) served as a positive control. The three negative controls are X 
DNA (S), human genomic DNA (Q, and 6X1 74 DNA (D). 



Karyotype Depiction of the Array Hybridization 
Pattern 

The inserts contained in the arrayed X clones 
have been mapped physically (Riles et al. 1993). 
The clones are arrayed in a random but known 
order on the array. Therefore, using the identity 
of each clone along with its physical map infor- 
mation, the pattern of hybridization to the yeast 
array can be represented in the form of a karyo- 
type of the yeast genome, as shown in Figure 2. 
The color of any segment of the ideogram repre- 
senting an individual chromosome on the karyo- 
type is directly determined by the ratio of red and 
green hybridization signals at the array positions 
of the corresponding clones. The lengths of the 
discrete colored segments of each chromosome 
correspond to the physical lengths of the yeast 



inserts. The chromosome seg- 
ments colored black represent ei- 
ther intervals of the genome that 
are not represented by clones in 
the library (90%> or false-negative 
hybridization signals on the array 
(10% i. Most of these false nega- 
tives are attributable to failures of 
the PCR amplification of the X 
clones, though occasional failures 
of the arraying process or nonuni- 
form surface preparation could ac- 
count for a small fraction of the 
false-negative signals. The large 
gap on chromosome 12 is the re- 
gion coding for ribosomal DNA 
that was not represented among 
the arraved clones. Genomic inter- 
vals represented, by overlapping 
clones were assigned a color based 
on the hybridization signals of 
only one of the overlapping 
clones, chosen at random. 

Note that in this representa- 
tion of a yeast karyotype, the larg- 
est six chromosomes are mainlv 
colored red. This indicates that 
most of the arraved clones that 
were mapped previously to these 
six large chromosomes hybridized 
primarily to the lissamine-labeled 
probe prepared from the corre- 
sponding purified chromosomes. 
Conversely, the smallest 10 chro- 
mosomes are mainly colored green 
in this image, matching the origi- 
nal CHEF gel isolation of the chro- 
mosomes used as the hybridization probe. The 
experiment was repeated with the yeast genome 
split into six discrete chromosome pools contain- 
ing 2-4 chromosomes per pool using CHEF gel 
electrophoresis. The chromosomes in each pool 
were extracted from the gel, amplified, and fluo- 
rescently labeled. The six chromosome pools 
were hybridized to six separate yeast arrays. 
Forty-four X clones gave a positive hybridization 
signal on all six arrays indicating that they con- 
tain yeast repetitive sequences (data not shown). 
These 44 clones and 10 clones with very weak 
hybridization signals were not included in the 
data set used to produce this karyotype. 

There were -40 anomalous clones, which ap- 
pear in this karyotype representation as green 
bands on the otherwise red chromosomes or red 
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Figure 2 Computer-generated ideogram repre- 
senting a karyotype of 5. cerevisiae, based on the 
normalized hybridization signals from the array 
shown in Fig. 1 . Note that the 6 largest chromo- 
somes are mainly red and the 10 smallest chromo- 
somes are mainly green. Black stripes represent in- 
tervals not represented by clones in the array or for 
which the corresponding clones gave false-negative 
hybridization signals. 

bands on the otherwise green chromosomes. 
Four randomly chosen examples of these anoma- 
lous clones were analyzed by hybridizing the 
clones to vertical strips cut from a Southern blot 
of CHEF gel-separated yeast chromosomes. In 
each case, the hybridization patterns of the 
anomalous clones corroborated the chromo- 
somal locations assigned by the microarray hy- 
bridization results (data not shown). Two clones 
that were thought to map to the 10 smallest chro- 
mosomes were found to hybridize preferentially 
to the probe representing the 6 largest chromo- 
somes and thus appear as anomalous red bands 
on the karyotype. Both hybridized to one of the 
six largest chromosomes on the Southern blot. 
Similarly, two clones that appear as anomalous 
green bands on the karyotype were found to hy- 
bridize to one of the 10 smallest chromosomes on 
the Southern blot. Thus, the anomalous clones 
are probably the result of sample tracking errors 
or, possibly, of errors in the published restriction- 
digest-based physical map on which the karyo- 
type representation was based (Riles et al. 1993). 

DISCUSSION 

The DNA microarray hybridization system re- 
ported here is conceptually and functionally 



similar to fluorescent in situ hybridization (FISH) 
to metaphase chromosomes, with three impor- 
tant differences. First, the target elements of the 
microarrays can, in principle, be any length or 
composition, from megabase YAC clones or mi- 
crodissected chromosome bands to individual 
cDNA clones, to short oligonucleotides. This ver- 
satility allows the user to choose characteristics, 
such as the mapping resolution and genetic com- 
plexity of each array element, to suit a particular 
application. Second, the hybridization signals are 
localized to discrete elements of known size and 
location, making them easier to identify and 
quantitate than the hybridization signals from 
irregularly shaped metaphase spreads. Third, mi- 
croarrays are more consistent and potentially 
amenable to automated production, hybridiza- 
tion, and data analysis than metaphase spreads. 

Arrays of DNA samples on porous mem- 
branes, for example, dot blots, have long been 
used as a basic tool in molecular biology. Dot- 
blot membranes are usually at least 8 x 12 cm in 
size, require the use of milliliter volumes of hy- 
bridization solution, and are limited, owing to 
autofluorescence and scattering, to radioactive, 
chemiluminescent, and colorimetric hybridiza- 
tion detection methods (Ross et al. 1992). Micro- 
arrays made on glass surfaces, on the other hand, 
can be mass-produced and are comparatively in- 
expensive, convenient, and compatible with 
fluorescent hybridization detection methods. 
Furthermore, a glass surface, when appropriately 
treated, has very low nonspecific binding of la- 
beled hybridization probes, resulting in lower 
backgrounds than are encountered typically with 
porous membranes. For hybridizations with very 
complex probes, the concentration of the labeled 
probe DNA is a limiting factor in the sensitivity 
of the assay. Minimizing the volume of the probe 
solution in a hybridization, by restricting the tar- 
get to a small area and by using a nonporous 
substrate, makes it practical to achieve very high 
probe concentrations. 

One important advantage of fluorescently la- 
beled probes is that, unlike most radioactive and 
chemiluminescent signals, fluorescent signals do 
not disperse and therefore allow for very dense 
array spacing. A unique, and probably the most 
important, advantage of fluorescent probes is 
that the hybridization signals from two or more 
differently labeled probes hybridized to the same 
target element can be detected separately. In this 
way, two-color hybridization detection allows for 
a direct and quantitative comparison of the 
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abundance of specific sequences between two 
probe mixtures that are hybridized competitively 
to a single anay. The absolute intensity of a hy- 
bridization signal at a particular element in an 
array can vary owing to experimental factors 
such as variations in the amount of DNA depos- 
ited on the array, variations in the hybridization 
or wash conditions between experiments, or 
variations in the hybridization characteristics of 
the different DNA sequences on the array. The 
ratio of the two signals at any element in an ar- 
ray, however, is relatively insensitive to these 
confounding factors because they affect both 
probe mixtures equivalently. This ratio therefore 
accurately reflects the relative abundance of the 
cognate sequence in the two probe samples. This 
is the principle underlying the technique of com- 
parative genomic hybridization (CGH), which is 
used to detect changes in the copy number of 
specific chromosomes or chromosomal regions 
(Kallioniemi et al. 1992). CGH is based on mea- 
suring the relative fluorescent hybridization in- 
tensities of two genomic-complexity hybridiza- 
tion probes, for example, probes representing ge- 
nomic DNA from normal and affected tissue 
samples, which are labeled with two distinct fluo- 
rophores and hybridized simultaneously to a 
metaphase spread. DNA microarray representa- 
tions of the human genome may provide a more 
convenient and higher resolution alternative to 
metaphase chromosomes for CGH. 

Cross-hybridization between related se- 
quences is an important problem faced by any 
hybridization-based assay, including the DNA 
microarray assay described here. Studies are now 
in progress to quantitate the extent of cross- 
hybridization between related sequences of vary- 
ing homology and length, in DNA microarray 
hybridizations. The stringency of hybridization 
and washing can be controlled by varying the salt 
concentration and temperature as in conven- 
tional membrane-based hybridizations. Cross- 
hybridization caused by repetitive sequences can 
be minimized by prehybridization of the probe or 
array with vast excess of unlabeled copies of the 
repetitive sequences. 

Alternative methods have been described for 
making microarrays of very short DNA se- 
quences, involving photolithography (Pease et 
al. 1994) or physical masking (Maskos and South- 
ern 1992) methods. These in situ synthesis meth- 
ods are inherently limited to low complexity ar- 
ray elements consisting of oligonucleotides. For 
complex-probe hybridizations, the specificity of 



hybridization is improved by using DNA frag- 
ments substantially longer than oligonucleo- 
tides. Moreover, the in situ synthesis approaches 
to array fabrication depend on prior knowledge 
of the sequence to be recognized by each array 
element. The approach described here makes mi- 
croarrays by transferring tiny volumes of DNA 
samples from microwell storage plates to a solid 
substrate. Thus, nucleic acids (or other mol- 
ecules) of virtually any length or any origin can 
be arrayed, and knowledge of their sequences is 
not required. 

The arrays used in these experiments do not 
represent the maximal achievable density of ele- 
ments. We have found that the spacing between 
the spots can be decreased by shrinking the con- 
tact area of the printing tip and by increasing the 
hydrophobicity of the glass surface. Microarrays 
with 100-jim feature size have been tested suc- 
cessfully in pilot experiments (data not shown). 
Assuming the projected availability of the appro- 
priate physically mapped human genomic clones 
(Hudson et al. 1995), arrays at 100-p.m spacing 
would allow for 10,000 discrete intervals of the 
human genome to be represented in a 1-cm 2 ar- 
ray. Such an array could be used for mapping at a 
resolution of <0.5 Mb. Experiments are in 
progress to explore the feasibility of such arrays. 

Our initial motivation for developing these 
microarrays arose from the need for abundant 
and inexpensive genomic arrays for genomic 
mismatch scanning (GMS) (Nelson et al. 1993), a 
method of genetic linkage analysis based on 
identification of the regions of "identity by de- 
scent" between affected relative pairs using a 
single complex-probe hybridization to an array 
of genomic clones. Experiments using these ar- 
rays to map quantitative trait loci in yeast by . 
GMS are currently in progress 0- deRisi, D. Lash- 
kari, L. Penland, L. McAllister, J. McCusker, R. 
Davis, and P.O. Brown, unpubl.). 

Microarrays of cDNA clones, prepared using 
the system described here, have been used for 
quantitative monitoring of gene expression pat- 
terns in Arabidopsis (Schena et al. 1995), S. cercvi- 
siae (D. Lashkari, J. deRisi, L. Penland, P.O. 
Brown, and R. Davis, unpubl.), and human tis- 
sues 0- deRisi, M. Bittner, P. Meltzer, L Penland, 
J. Trent, and and P.O. Brown, unpubl.). We an- 
ticipate that DNA microarrays of the kind de- 
scribed here will be useful in additional applica- 
tions for which conventional dot blots, high- 
density gridded arrays on porous membranes, or 
FISH are currently used. These potential applica- 
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tions include c mparative genomic hybridiza- 
tion (Kallioniemi et al. 1992), sequencing by hy- 
bridization (Drmanac et al. 1993), physical map- 
ping of cloned or amplified sequences (Billings et 
al. 1991), and economical distribution of re- 
agents for integrated genetic and physical map- 
ping based on a common set of arrayed clones 
(Zehetner and Lehrach 1994). 

METHODS 

Amplification of Target DNA Elements 

The anay elements were prepared from physically mapped 
x cl nes (Riles et al. 1993). The k clones were amplified 
using randomly primed polymerase chain reaction (PCR) 
based n published and unpublished protocols (Bohlandei 
et al. 1992; S. Nelson, unpubU. The phage lysates were 
amplified in a 10-m-I PCR reaction using 5 yM final concen- 
trate n of primer A (GCTATCTTCAAGATCANKNNNN), 
200 *lm dNTPs. and 1 unit of Taq polymerase. Round A 
consisted of five cycles at 94*C for 1 min. 25 B C for 1.5 min, 
25-72*C over 7 min. and 72*C for 3 min using Taq poly, 
merase (BMB}. For round B, the reaction volume was 
brought up to 100 mJ for a final concentration of 2 »im of 
primer B (GCTATCTTCAAGATCAj, 200 n* dNTPs, and 4 
units f Taq polymerase. Round B consisted of 30 cycles of 
94*C f r 1 min, S6*C for 2 min, and 7TC for 3 min. The 
amplification was performed in 96-well plates using crude 
phage lysates as the templates, resulting in an amplified- 
ti n f both the 35-kb X vector and the 5-kb to 15-kb veast 

m 

insen sequences as a distribution of PCR products between 
250 bp and 1500 bp in length. 

The PCR products were purified and transferred into 
TE (10 mM Tris, 1 mM EDTA at pH 8.0) buffer using Sepha- 
dex G50 gel filtration (Pharmacia) and evaporated to dry- 
ness at room temperature overnight. Each of the 864 am- 
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plified X clones was rehydrated in 15 *1 of 3x SSC 
(20 x 55C « 3 m NaCl. 0.3 u Na, citrate) in preparation foi 
spomng onto the glass under n rmal room temperature 
conditions. 



Preparation of DNA Microarrays 

The microarrays were fabricated on poly-i-lysine coated 
microscope slides (Sigma). A custom-built arraying ma- 
chine, consisting of four tweezer-like printing tips 
mounted 9 mm apart on a computer-controlled robotic 
stage (Shalon 1996), loaded 1 m>1 of the concentrated PCR 
product directly from corresponding clusters of four wells 
of 96-well storage plates and deposited -5 nl of each 
sample onto each of 40 slides. Surface tension loaded the 
sample into the printing tip directly from the miaowell 
plate and held the sample in the tip during the printing 
operation. Printing was achieved by lightly tapping the tip 
against the glass surface. The open-capillary design al- 
lowed for rapid nnsing and drying of the tips between 
samples. Figure 3 shows the layout of the arraying ma- 
chine. Figure 4 shows a detailed view of the four printing 
tips and the staggered printing partem on the microscope 
slides. Adjacent samples were spotted 380 urn apart on the 
slides. After each set of four samples was printed onto 40 
slides, the printing tips were rinsed with a tet of water for 
2 sec and then dried by lowering the tips onto a sponge for 
2 sec. The process was repeated for all 864 samples and 
eight control spots. 

After the spotting operation was complete, the slides 
were rehydrated in a humid chamber at room temperature 
for 2 hr, baked in an 80*C vacuum oven for 2 hr, then 
rinsed in 0.1% sodium dodecyl sulfate (SDS) to remove 
unadsorbed DNA. To reduce nonspecific adsorption of the 
labeled hybridization probe to the poly-L-lysine coated 
glass surface, the slides were treated with succinic anhy- 
dride. One gram of succinic anhydride was dissolved in 
100 ml of l-methyl.2-p>Trolidinone and then 100 ml of 

0.2 m boric acid (pH 8.0) was 
added. The arrays were soaked in 
this solution for 10 min and then 
rinsed in distilled water four 
times for 5 min each. Immedi- 
ately before use. the arrayed DNA 
elements were denatured by plac- 
ing the slide in distilled water at 
90*C for 2 min. 
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Flgur 3 The layout of the arraying machine. All motions 
control. For more details f the arraying machine, see 
cmgm.stanford.edu/pbrown. 
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Amplification and Labeling 
of Hybridization Probe 

The 16 chromosomes of Saccharo- 
myces cerrvisiae were separated us- 
ing a contour-clamped homoge- 
neous electric field (CHEF) aga- 
r se gel apparatus (Bio-Rad) (Qui 
et al. 1986). The 6 largest chromo- 
s mes were isolated in one gel 
slice and the smallest ten chro- 
mosomes in a second gel slice. 
The DNA fr m each slice was re- 
covered using a gel extraction kit 
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Microscope slide 

Figure 4 A close-up view of the four open- 
capillary printing tips. The tips are 9 mm apart and 
fit into four adjacent wells of a standard microwell 
plate and print arrays in a staggered fashion on mi- 
croscope slides. For more details of the printing tips, 
see web page http://cmgm.stanford.edu/pbrown. 



(Qiagen) and randomly amplified in a manner similar to 
that used in amplifying the target k clones (Grothues et al. 
1993). The main difference between this amplification 
procedure and the one used for the k array elements is a 
filtration step between rounds A and B to remove primer- 
dimers and the use of a random 9-mer 3' end on primer A. 
Following amplification. 2.5 M>g of each of the amplified 
chromosome pools were separately random-primer labeled 
using Klenow polymerase (Amersham) with a lissamine- 
conjugated nucleotide analog (DuPont NEN) for the pool 
containing the 6 largest chromosomes and with a fluores- 
cein-con jugated nucleotide analog f BMB) for the pool con- 
taining the smallest 10 chromosomes. The two fluores- 
cent-labeled pools were mixed and concentrated using an 
ultrafiltration device (Amicon). 



Hybridization 

Five micrograms of the hybridization probe, consisting of 
both chromosome pools in 7.5 of TE, was denatured in 
a boiling water bath and then snap-cooled on ice. Concen- 
trated hybridization solution (2.5 was added to a final 
concentration of 5x SSC/0.1% SDS. The entire 10 mJ of 
probe solution was transferred to the array surface, covered 
with a coversiip. placed in a custom-built single-slide hu- 
miditv chamber, and incubated in a 60*C water bath for 12 
hr. The custom-built waterproof slide chamber has a cavity 
just slightly bigger than a microscope slide and was kept at 
100% humidity internally by the addition of 2 wJ of water 
in a comer of the chamber. The slide was rinsed in 5 x 
SSC/0.1% SDS for 5 min and then in 0.2x SSC/0.1% SDS 
for 5 min. All rinses were at room temperature. The array 
was then air dried, and a drop of antifade (Molecular 
Probes) was applied to the array under a 24-mm x 30-mm 
coversiip in preparation for scanning. 



Detection and Analysis 

A cust m-built laser scanner was used to detect the two 



c lor flu rescence hybridization signals from 1.8- 
cm x 1.8-cm arrays at 20-iun resolution. The glass sub- 
strate slide was mounted on a c mputer-controlled. two. 
axis translation stage (PM-500. Newport. Irvine. CA) that 
scanned the array over an upward-facing microscope ob- 
jective (20 x, 0.75NA Fluor. Nikon. Melville. NY) in a bi- 
directional raster pattern. A water-cooled Argon/Krypton 
laser (Innova 70 Spectrum. Coherent. Palo Aito. CA). op- 
erated in multiline mode, allowed for simultaneous speci- 
men illumination at 488.0 nm and 568.2 nm. These two 
lines were isolated by a 488/568 dual-band excitation filter 
(Chroma Technology, Brattleboro. VT). An epifluores- 
cence configuration with a dual-band 488/568 primary 
beam splitter (Chroma) excited both fluorophores simul- 
taneously and directed fluorescence emissions toward the 
two-channel detector. Emissions were split by a secondary 
dichroic mirror with a 565 transition wavelength onto two 
multialkali cathode photomultiplier tubes (PMT; R928, 
Hamamatsu, Bridgewater. NJ), one with an HQS3S/50 
bandpass barrier filter and the other with a D630/60 band- 
pass barrier filter (Chroma). Preamplified PMT signals were 
read into a personal computer using a 12-bit analog-to- 
digital conversion board (KT1-834, Analog Devices, Nor- 
wood, MA), displayed in a graphics window, and stored to 
disk for further rendering and analysis. The back aperture 
of the 20 x objective was deliberately underfilled by the 
illuminating laser beam to produce a large-diameter illu- 
minating spot at the specimen (5-um to 10-p.m half- 
width). Stage scanning velocity was 100 mm/sec. and PMT 
signals were digitized at 100 psec intervals. Two successive 
readings were summed for each pixel, such that pixel spac- 
ing in the final image was 20 tun. Beam power at the 
specimen was -5 mW for each of the two lines. 

The scanned image was despeckJed using a graphics 
program (Hijaak Graphics Suite) and then analyzed using 
a custom image gridding program that created a spread- 
sheet of the average red and green hybridization intensi- 
ties for each spot. The red and green hybridization inten- 
sities were corrected for optical cross talk between the fluo- 
rescein and lissamine channels, using experimentally 
determined coefficients. 
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